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Universal properties of three-dimensional random-exchange quantnm antiferromagnets 

D.-R. Tan^ and F.-J. Jiang^ 0 

^Department of Physies, National Taiwan Normal University, 88, Sec.4, Ting-Chou Rd., Taipei 116, Taiwan 

The thermal and ground state properties of a class of three-dimensional (3D) random-exchange 
spin-1/2 antiferromagnets are studied using hrst principles quantum Monte Carlo method. Our 
motivation is to examine whether the newly discovered universal properties, which connect the 
Neel temperature and the staggered magnetization density, for the clean 3D quantum dimerized 
Heisenberg models remain valid for the random-exchange models considered here. Remarkably, 
similar to the clean systems, our Monte Carlo results indicate that these universal relations also 
emerge for the considered models with the introduced antiferromagnetic randomness. The scope 
of the validity of these universal properties for the 3D quantum antiferromagnets is investigated as 
well. 

PACS numbers: 


I. INTRODUCTION 

Heisenberg-type models have been studied extensively 
using both analytic and numerical methods during the 
last two decades [ll-[T^. This is mainly due to the fact 
that these models can qualitatively or quantitatively de¬ 
scribe experimental data. In real materials, impurities of¬ 
ten play important roles in their properties. As a result, 
generalized Heisenberg-type models such as taking into 
account the effects of (non)magnetic impurities have been 
explored in great detail as well fisl - f^ . In additional to 
impurities, quench disorder effects, namely randomness 
in the strength of the antiferromagnetic coupling, also 
attract a lot of theoretical interest [28l - l^ . Indeed, in¬ 
vestigation associated with impurity and disorder effects 
for antiferromagnets has led to several exotic (theoreti¬ 
cal) results. Two notable such examples are the anoma¬ 
lous Curie constant and violation of Harris criterion 
(3^ . While most of these studied are related to two- 
dimensional models, experimental results, like those of 
TlCuCH . have triggered theoretical attention to 

higher dimensional systems as well (stI - I^ . It should be 
pointed out that for three spatial dimensions which is the 
upper critical dimension, the relevant universal physical 
quantities close to (second order) quantum phase tran¬ 
sitions, namely the critical exponents are described by 
the corresponding mean-field values (Beside the leading 
exponents, there are logarithmic corrections to the scal¬ 
ing as well). Still, investigating universal relations that 
do not depend on the microscopic details of the systems 
is an interesting research topic. Such studies are also 
practically useful from a experimental point of view. 

Recently, close to quantum critical points (QCPs), sev¬ 
eral universal relations between the Neel temperature 
Uv and the staggered magnetization density Mg are es¬ 
tablished for three-dimensional (3D) spin-1/2 dimerized 
antiferromagnets both theoretically and experimentally 
(iO - li^l . For instance, as functions of Mg, the physical 
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quantities Tpf/J and T/v/T* show universal dependence 
on Mg [dlj. Here J and T* are the averaged strength of 
antiferromagnetic couplings and the temperature where 
peak of the uniform susceptibility as a function of tem¬ 
peratures takes place, respectively. Notice these univer¬ 
sal relations simplify to linear ones in the vicinity of 
quantum critical points. Interestingly, for 3D random- 
exchange systems, it is demonstrated in Refs. [13, 
that the long-range antiferromagnetism is robust against 
both the box-like and the singular disorder distributions. 
In addition, a linear dependence of T/v/J on Mg close to 
the correspondii^data points of the clean systems is ob¬ 
served in Refs. [31^ as well (In this study observables 
with a overline on them stand for the disordered average). 
This is quantitatively different from the scenario of clean 
3D dimerized Heisenberg models. Notice for the random- 
exchange models considered in Refs. [H, E^, the related 
quantum phase transitions do not occur at finite random¬ 
ness in the relevant parameter spaces. Hence it would be 
extremely interesting to investigate whether the univer¬ 
sal relations between T/v and Mg, which are found for 
the clean 3D quantum dimerized systems, remain valid 
for 3D spin-1/2 random-exchange models with the cor¬ 
responding quantum critical points taking finite values 
in the related parameter spaces. Motivated by this, here 
we investigate a class of 3D quantum random-exchange 
Heisenberg models which undergo quantum phase transi¬ 
tions with randomness of finite magnitude. These models 
will be introduced explicitly later. Remarkably, the uni¬ 
versal relations observed in Ref. [4l| are valid for the 
models considered in this study. Specifically, for our 
models, both the quantities Tjv/J and T^/T* do show 
universal dependence on Mg. We have further studied 
the scope of the validity of these universal relations and 
concluded that such properties are justified within each 
individual category (explained later). 

The rest of this paper is organized as follows. After the 
introduction, the studied models, the employed technical 
methods, the used randomness, as well as the observables 
are described in section 2. Section 3 contains the detailed 
numerical results and discussion. Finally we concludes 
our investigation in section 4. 
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II. MICROSCOPIC MODEL, BOND 
RANDOMNESS AND OBSERVABLES 


The studied 3D random-exchange quantum Heisenberg 
models are given by the Hamilton operators 

h = Y, S. -Sj + Y, 4 • 4 , ( 1 ) 

(b) d'i') 


where Jij and J[iy are the antiferromagnetic exchange 
coupling connecting nearest neighbor spins {ij) and {i'j'), 
respectively, and Si is the spin-1/2 operator at site i. 
For simplicity, here we will call Jij and the (an¬ 
tiferromagnetic) bonds. Fig. [1] demonstrates the mod¬ 
els described by Eq. 0 in the pictorial form. Notice 
in this study the models shown in the top and bot¬ 
tom panels of Fig. [T] will be called the (random) ladder- 
and staggered-dimer models, respecitvely. The random¬ 
ness to the antiferromagnetic coupling strength is intro¬ 
duced as follows. First of all, the strength of the anti¬ 
ferromagnetic couplings Jij are set to 1.0 for both mod¬ 
els. Second, J[, j, takes the value of Jc (1.0 -I- D) and 
Jc( 1.0 — D) randomly with probabilities P and 1 — P, 
respectively. Here 0.0 < D < 1.0 and 0.0 < P < 1.0. 
The Jc appearing above is the critical coupling of the 
clean system which is given by Jc = 4.0128 {Jc = 4.6083) 
for the random ladder-dimer (staggered-dimer) quantum 
Heisenberg model [3^. A similar (2D) model has been 
considered in as well. Based on the method of in¬ 
ducing randomness, there exists a critical Pc at which 
a quantum phase transition from the antiferromagnetic 
long-range order to a disordered phase occurs. Here we 
will focus on the case of D = 0.5. To determine the 
Neel temperature Tjv and the staggered magnetization 
density Mg for the considered models with the employed 
random ness, the observables staggered structure factor 
5 '(7r, tt), spin stiffness and second Binder ratio Q 2 are 
calculated in our simulations. The quantity iS'(7r,7r) is 
dehned by 

w:^ = (2) 

on a finite cubical lattice with linear size L. Here = 
with Sf being the third-component 
of the spin-1/2 operator Si at site i. In addition, ^ is 
given as 

(3) 


where Wi is the winding number in i direction and {3 is 
the inverse temperature. Finally, the observable Q 2 takes 
the form 


Q 2 


{{mlY) ■ 


( 4 ) 


With these observables, T/v and Mg can be determined 
with high precision. 
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FIG. 1: The three dimensional random ladder- (top panel) 
and staggered-dimer (bottom panel) quantum Heisenberg 
models considered in this study. 


III. THE NUMERICAL RESULTS 


To understand whether universal relations between Tjv 
and Mg, similar to those of clean dimerized systems, ap¬ 
pear for the studied 3D quantum Heisenberg models with 
the introduced bond randomness, we have carried out 
a large-scale Monte Carlo simulation using the stochas¬ 
tic series expansion (SSE) algorithm with very efficient 
loop-operator update 14711 . In particular, the /3-doubling 
scheme described in |48l | is used in our investigation in 
order to access the ground state properties in an efficient 
way. Eurthermore, each random configuration is gener¬ 
ated by its own seed and several hundred (for obtaining 
ground state properties) to several thousand (for the de¬ 
termination of thermal properties) randomness realiza¬ 
tions are produced. The potential systematic uncertain¬ 
ties due to thermalization, Monte Carlo sweeps within 
each randomness realization, as well as the number of 
configurations used for disordered average are examined 
by carrying out many trial simulations. The resulting 
conclusions based on these trial simulations are consis¬ 
tent with those presented in this study. 


A. The determination of Mg 


The reach of the ground state values of S{-k, tt) for sev¬ 
eral considered P and L is shown in fig. [2] for both the 
random ladder- and staggered-dimer models. Further¬ 
more, the determination of Mg is done by extrapolating 
the related finite volume staggered structure factors to 
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FIG. 2: Convergence of ^(TrjTr) to the ground state values, 
for several considered P and various box sizes L, for the ran¬ 
dom ladder- (top panel) and staggered-dimer (bottom panel) 
models. The solid lines are added to guide the eye. 


FIG. 3 : L dependence of the staggered structure factors 
^(TrjTr), at several considered values of P, for both the ran¬ 
dom ladder- (top panel) and staggered-dimer (bottom panel) 
models. The dashed lines are added to guide the eye. 


the corresponding bulk results, using a polynomial fit of 
the form a^+ai/L + a 2 /L^+a ^/(The bulk Mg is given 
by y/oo). For both the studied models, the L-dependence 
of S'(7r,7r) for several considered P is depicted in fig. El 
The obtained extrapolating values of Mg are shown in 

fig. El 


B. The determination of Tn 

The employed observables for calculating Tjv are 
{]h/J)L as well as < 52 . Notice a constraint standard 
finite-size scaling ansatz of the form (I -|- 6 oT““)( 6 i -I- 
+ + (or bi + b 2 tL^/^ + b 3 itL^/''f + ... 

in some cases), up to fourth order in , is used to fit 
the data. Here 6 ^ for f = 0,1, 2,... are some constants and 
f _ t-Tm 1300 J 1 demonstrated in Ref. H, [S that 

such an analysis procedure leads to accurate determina¬ 
tion of Tm- The data of (JT^/ J)L for the random ladder- 
dimer model with P = 0.3 (top panel) and the random 
staggered-dimer model with P = 0.05 (bottom panel) are 
shown in fig. [5] The values of T/v obtained from J)L 
are listed in table 1. Furthermore, the quoted errors ap¬ 
pearing in table 1 are conservatively estimated, based 
on the standard deviations obtained from a bootstrap 
method used for the fits. Finally, applying a similar anal¬ 
ysis to Q 2 leads to consistent results of Taj with those in 
table 1. Data points of Q 2 for the random ladder-dimer 
model with P = 0.5 and the random staggered-dimer 
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FIG. 4: Ms as functions of the considered values of P for 
both the studied models. While the top panel is for the ran¬ 
dom ladder-dimer model, the results of Ms for the random 
staggered-dimer model is depicted in the bottom panel. The 
dashed lines are added to guide the eye. 
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FIG. 5: {'pl/J)L as functions of T/J for L —12, 16, 20, 24, 
28, 32, and 36 for both the studied models. J is 1.0 in our 
calculations. While the top panel is for the random ladder- 
dimer model with P = 0.3, the results of {'p^/J)L for the 
random staggered-dimer model with P — 0.05 is depicted in 
the bottom panel. The dashed lines are added to guide the 
eye. 




T/J 


FIG. 6: Q 2 as functions ol T/J for L =12, 16, 20, 24, 28, 
32, and 36 for both the studied models. J is 1.0 in our 
calculations. While the top panel is for the random ladder- 
dimer model with P = 0.5, the results of Q 2 for the random 
staggered-dimer model with P — 0.2 is depicted in the bottom 
panel. The dashed lines are added to guide the eye. 


model with P = 0.2 are presented in fig. |6l In additional 
to the method of finite-size scaling, T/v can also be cal¬ 
culated by studying the singular behaviour of Mg when 
approaching Tn from T < T]\[. We have performed such 
investigation for the random ladder-dimer (staggered- 
dimer) model with P = 0.05,0.55 {P = 0.1,0.4). The 
obtained results of Tjv are in remarkable agreement with 
those shown in table 1. See figs. [7] and |8] for the details. 


C. The determination of T* 

In additional to the universal relation between T]\j/J 
and Mg, it is demonstrated in Ref. [4l| that there exists a 
universal connection between the quantities T^/T* and 
Mg as well. Here T* is the temperature where peak of 
the uniform susceptibility (x„ = (^(E,-5f)2)),asa 
function of the temperatures, takes place. For both the 
studied models, the values of T* for the considered P are 
estimated from the related data of L = 12 and L = 24. 
The numerical values of Xu as functions of T/J for some 
values of P and L = 12 are shown in fig. [H 


p 

Tn/J 

P 

Tn/j 

0.025 

1.01030(48) 

0.05 

1.02030(53) 

0.05 

0.99909(52) 

0.1 

0.99308(48) 

0.1 

0.97550(40) 

0.2 

0.93097(49) 

0.2 

0.92208(40) 

0.3 

0.85645(44) 

0.3 

0.85703(50) 

0.4 

0.76297(53) 

0.35 

0.81891(39) 

0.45 

0.70586(66) 

0.45 

0.72636(39) 

0.5 

0.63896(56) 

0.5 

0.66832(44) 

0.55 

0.55730(21) 

0.55 

0.59868(62) 

0.6 

0.44937(31) 

0.6 

0.51146(32) 



0.625 

0.45649(28) 




TABLE I: Numerical values of Tn / J determined from the ob¬ 
servable (j}^/J)L for various P. J is 1.0 in our calculations. 
The first two columns and the last two columns are associ¬ 
ated with the random ladder- and staggered-dimer models, 
respectively. 




















5 




FIG. 7: Ms/y/i as functions oi T/J for the random ladder- 
dimer model with P = 0.05 (top panel) and P = 0.55 (bottom 
panel). J is 1.0 in our calculations. The data points are 
determined with a up to third order polynomial formula in 
1/L and the quoted errors are estimated directly from the 
fits. The solid lines are obtained using the results of the fits 
(leading singular behaviour ansatz). 




FIG. 8: Ms/ Vs as functions of T/ J for the random staggered- 
dimer model with P = 0.1 (top panel) and P = 0.4 (bottom 
panel). J is 1.0 in our calculations. The data points are 
determined with a up to third order polynomial formula in 
1/L and the quoted errors are estimated directly from the 
fits. The solid lines are obtained using the results of the fits 
(leading singular behaviour ansatz). 


D. The universal relations 

After having obtained Tjv, J, Mg, and T*, as a first 
step to examine whether universal relations, as those es¬ 
tablished for regular dimerized systems, would appear 
for the considered models with the employed random¬ 
ness, we study T/v as functions of Mg for both models. 
The results are shown in fig. [TUI As one can see in fig. (TUI 
no clear signal of a universal relation is found in the fig¬ 
ure. Remarkably, if the quantity Tjq/J is investigated 
as functions of Mg , then indeed the data points of both 
models fall on top of a universal curve. The result is 
demonstrated in fig. 1111 Similarly, a universal behaviour 
between Tn/T* and Mg is observed as well (fig. [T2ll . 


IV. DISCUSSIONS AND CONCLUSIONS 

Inspired by the universal relations between the Neel 
temperature Tv and the staggered magnetization density 
Mg for the clean 3D dimerized systems, in this work we 
study these universal behaviour of Tv and Mg for a class 



T/J 


FIG. 9: xL of the random ladder-dimer model as functions 
of T/J for some values of P and P = 12. J is 1.0 in our 
calculations. 


of 3D random-exchange quantum Heisenberg models. A 
notable finding here is that these universal properties 
are even valid for the considered models with the intro¬ 
duced randomness. Our results indicate that the scope of 
the validity of these universal properties for 3D quantum 
antiferroinagnets is very general. One interesting ques- 
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FIG. 10: Tn as functions of Ms for both the considered models 
with the introduced randomness. 
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dition, the data points associated with clean dimerized 
ladder model do not form a universal curve with those 
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FIG. 12: Universal behaviour of Tn/T* as functions of Ms 
for the considered models with the introduced randomness. 
The corresponding results of Tn are shown in the figure as 
well. The data points of Tn/T* contain (around) one percent 
uncertainties due to the errors of T*. 

data points presented in fig.[TTJ As a result, these univer¬ 
sal properties are valid within individual categories, such 
as the models in fig. [T] with D = 0.5 which are investi¬ 
gated in detail here. Nevertheless, it is remarkable that 
the universal relations originally observed for the clean 
models persist even for systems with randomness. 


FIG. 11: Universal behaviour of Tjv/J as functions of Ms for 

the considered models with the introduced randomness. Acknowledgments 


tion is to examine whether the data points obtained with 
other values of D will fall on top of the universal curves 
shown in figs. Illl and ll2l Our preliminary results related 
to Tm/J and Mg show that this is not the case. In ad- 
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